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NASA’s Solar Sail Demonstration project partner L’Garde is currently assembling a flight- 
like sail assembly for a series of ground demonstration tests beginning in 2015. For future 
missions of this sail that might validate solar sail technology, it is necessary to have an 
accurate sail thrust model. One of the primary requirements of a proposed potential 
technology validation mission will be to demonstrate solar sail thrust over a set time period, 
which for this project is nominally 30 days. This requirement would be met by comparing a 
L’Garde-developed trajectory simulation to the as-flown trajectory. The current sail 
simulation baseline for L’Garde is a Systems Tool Kit (STK) plug-in that includes a custom- 
designed model of the L’Garde sail. The STK simulation has been verified for a flat plate 
model by comparing it to the NASA-developed Solar Sail Spaceflight Simulation Software 
(S5). S5 matched STK with a high degree of accuracy and the results of the validation 
indicate that the L’Garde STK model is accurate enough to meet the potential future 
mission requirements. Additionally, since the L’Garde sail deviates considerably from a flat 
plate, a force model for a non-flat sail provided by L’Garde sail was also tested and 
compared to a flat plate model in S5. This result will be used in the future as a basis of 
comparison to the non-flat sail model being developed for STK. 


Nomenclature 

A = sail area 

B b = back non-Lambertian coefficient 1 9 

Bf = front non-Lambertian coefficient 1,9 

F x = force component along sail x body axis 

F z = force component along sail z body axis 

P = solar pressure at one astronomical unit (AU) 

f = reflection coefficient 1 

R p = propulsive reflectivity coefficient 9 

s = fraction of specular reflection 1 

s b = back surface emissivity 1 ' 9 

£f = front surface emissivity 1 ' 9 

= sun incidence or pitch angle 1 9 
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9 p = propulsive zenith 9 

p d = total reflectivity 9 

a 2 = simplified optical constant for 1 st order tensor 

a 3 = simplified optical constant for 2 nd order tensor 

f = sail to sun unit vector in sail body frame 

Fc = sail force for Generalized Sail Model 

Jij k = tensors for Generalized Sail Model where n is order and ijk is index 

I. Introduction 

A s future missions studying the solar system are designed based on scientific objectives in areas such as 
planetary science, heliophysics, and interplanetary exploration, solar sails provide the ability to reach some 
orbits and locations either difficult to achieve with conventional chemical propulsion, or at a reduced overall 
cost. Solar sails produce thrust solely through the use of Solar Radiation Pressure (SRP) from the Sun, giving them 
the ability to operate in space with continuous thrust for long periods of time and without the need for propellant. 
The fundamental process from which this thrust is derived is the momentum transfer of photons reflecting off the 
sail surface. This process has been studied extensively, especially since it disturbs the orbits of all spacecraft 
exposed to the sun. Many solar sail studies and some planned missions have resulted from the understanding of this 
readily available source of thrust. At NASA specifically, these include the proposed but ultimately scrapped 
Halley’s Comet rendezvous in the 1970s 1 2 and the use of Mariner 10’s solar panels as solar sails for attitude 
control. 3 More recently, NASA funded ground demonstrations by ATK and L’Garde of 20-meter, fully-integrated 
sail systems in 2005 4 ' 5 (as seen in figure 1), and also NanoSail-D, a 3U Cubesat that deployed a 10 square meter sail 
in 2011 to test the use of a small sail for Low Earth Orbit (LEO) end-of-life de-orbit. 6 Also, the MESSENGER 
mission, launched in 2004 and continuing operation today, used SRP for trajectory and attitude control in reaching 
orbit around Mercury. 7 The Japanese Aerospace Exploration Agency (JAXA) also has been involved in developing 
and demonstrating the use of solar sails, and in 2010 deployed and launched the solar sail thrust demonstration 
spacecraft 1KAROS to Venus. 8 Even with the benefits of using a propellant-less source of continuous in-space 
thrust, because of the difficulties involved in manufacturing solar sail thin films and structures, along with the risks 
in deploying such huge structures in space and also those inherent in using an untested technology, solar sails 
currently still reside at the “proof-of-concepf ’ level for main spacecraft propulsion. 



Figures la and lb. ATK (left) and L’Garde (right) ground demonstration sail systems in the Glenn Research 
Center Plum Brook space environment simulation chamber . 4 ' 5 

Following the success of NanoSail-D and the ATK and L’Garde ground sail system demonstrations, NASA is 
funding the Solar Sail Demonstration Project as part of the Technology Demonstration Mission (TDM) program 
with sail manufacturer L’Garde. For this project, a set of requirements was established to support future Missions. 
The requirements included that a future mission of the sail be able to perform “a 30-day maneuver sequence [...] 
with less than 20 percent error by comparison with the predicted final velocity and position relative to Earth.*” For a 
future mission, this requirement is planned to be met by comparing the as-flown trajectory to a simulated trajectory 
based on the projected performance of the sail. The simulation L’Garde will use for this purpose is an AGI Systems 
Tool Kit (STK) “plug-in” model of the sail. At the time of the writing of this paper, the STK plug-in consisted of a 


* Mission requirements from Solar Sail Demonstration Project peer review May 2013. 
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flat plate sail with the L’Garde optics model. In the future, L’Garde plans to improve the solar sail model to include 
the model of their sail that incorporates shape. In this paper, we validate the flat plate model of L’Garde by 
comparing it to with the NASA-developed Solar Sail Spaceflight Simulation Software (S5). 9 We also use the shaped 
sail model in S5 to model the shape of the L’Garde sail and compare it to a flat plate model for future reference. 

II. Technical Approach 

The comparison between the L’Garde STK simulation and the corresponding simulation in S5 encompasses a 
number of factors. Each software tool incorporates different options for gravity models, SRP models, integrators, 
propagators, coordinate systems, and attitude representation. 

A. L’Garde STK Simulation 

The solar sail trajectory was simulated by L’Garde using STK Astrogator with a plug-in developed by Applied 
Defense Solutions (ADS) containing the solar sail propagator. Included along with the plug-in are the STK high 
precision orbit propagator (HPOP v8.1.1), which uses a gravity model ofEarth to degree and order 21, along with a 
point mass model of the Sun, and a point mass model of the Moon. The SRP force model used was that developed 
by Derbes and Lichodziejewski 1 " through tests on prototype sail materials. This model describes the forces in terms 
of a body coordinate frame as shown in Figure 2. 


(+) "Sun Incidence" 
about Ysail 


Figure 2. Diagram showing sail body axes and sun incidence angle 

The force developed by the sail from the SRP is given in the transverse and normal directions in Eq. (1) and Eq. 
(2) respectively below, both of which are functions of sun incidence angle 0j. P is the solar pressure at 1 AU of 
4.511918 X 10 -6 Pa and the other constants are defined in the appendix. 

F x — PAcos (0j) sin(0j) — P p PAcos(0i)sin(0 p ) (1) 

F z = -PAcos 2 (0*) - R p PAcos(e{) cos(0 p ) - (1 - p d )PAcos(6 i ) ( fB f^ b ) (2) 

The attitude of the sail is defined as in Derbes, 11 where top, sun incidence, and flat spin angles are a 3-2-3 Euler 
angle sequence about the body axes as are shown in Figure 2, starting from an initial orientation with the z body axis 
pointed at the Sun and the x body axis parallel to the ecliptic plane. Also, the propagation in STK used the Runge- 
Kutta 7/8 (Fehlberg) method. 
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For this study, the mass of the L’Garde sailcraft was assumed to be 34.641 kg, resulting in a characteristic 
acceleration of 0.289 mm/sec A 2. 

B. S5 Simulation 

The simulation of the relevant part of the Sunjammer mission using S5 allowed a selection of different options. 
The gravity model available in S5 was the J2 model of Earth, along with point mass gravity models of the Sun and 
the Moon. The SRJP models available in S5 are both the Mclnnes model 1 and the Derbes and Lichodziejewski 
model as above. In terms of similar body axes as in Figure 2, Eq. (3) and (4) give the Mclnnes model force in the 
transverse and normal directions. S5 uses a value of SRJP at 1 AU of 4.571176 X 10 -6 Pa in the Mclnnes SRP 
model, the forces are again functions of sun incidence angle 0 ; . In order for the Mclnnes model to match the Derbes 
model, we added a gain constant k to the 0, in the sine term of Equation (3) to make it match 0 p in the second term of 
Equation (1) (see Appendix for value of k). In the Derbes SRP model in S5 the value of SRP at 1 AU is consistent 
with that used in STK. 

F x = —PA(1 — rs ) cos(0j) sin(/c0j) (3) 

F z = — Pi4(l — rs)cos 2 (8i ) — PAB f (l — s)f cosCfy) — PA(1 — f)( E/g/ EbBb )cos (0j) (4) 

‘ £ f+ £ b 

The sail attitude was defined to match that in the L’Garde STK plugin, and the propagation in S5 used the Runge- 
Kutta 7/8 (Fehlberg) method also. 

C. Comparison 

After propagating from the same initial state using their respective force models, the final state at the end of the 
30 day maneuver sequence was output from both STK and S5 for comparison. Values of percent error were then 
calculated from the two state vectors to be consistent with the requirement as stated. Two different methods of 
calculating percent error were used for comparison as a final decision on the method of percent error calculation for 
the mission requirement (less than twenty percent error on position and velocity) has not yet been made. The first 
method uses vectors in the ICRF 12 coordinate frame and calculates percent error as in Eq. (5) and Eq. (6). 

Position Percent Error = 100 * H r ^2 t stk 2 II (5) 

\\ r STK2\\ 

Velocity Percent Error — 100 * 11^2 (6) 

\\ v STK2\\ 

Here r S s 2 * s l ^ le vector to the final position of the sail calculated using S5, r ST K 2 * s th e vector to the final position 
of the sail calculated using the STK plug-in, v S52 is the velocity of the sail at the end of the maneuver sequence 
calculated using S5, and finally, v ST k 2 is the velocity of the sail at the end of the maneuver calculated using STK. 

The second method for calculating the percent error between the two trajectories is again in the ICRF frame; 
however, the error values are instead calculated as in Eq. (7) and Eq. (8). 

Position Percent Error = 100 * l|Ar ,f 5 „ Ar5 J yl1 (7) 

IIArsTidl 

Velocity Percent Error = 100 * l|Al? ,f 5 _ (8) 

IIA^stwII 

In these two equations, A r S5 = r S5 2 — r S5l and A r STK = r s - n<2 — r STKl are the differences between initial and 
final position in S5 and STK respectively, and Av S5 = v S5 2 — v S5l and Av STK = v S tk 2 ~ Vstk 1 are the differences 
in initial and final velocities. Also, note that r S5l = r STKl and 1 = v ST k v as both sail simulations start from the 
same initial state. 

These two error methods were also used to estimate the effects of attitude error on the ability to achieve the 
mission requirement of 20% error on position and velocity. For this estimation the first two angles in the 3-2-3 Euler 
angle sequence (top spin about the z body axis and sun incidence about the y body axis respectively) were varied 
±2° in the S5 simulation. This variance is from the trim conditions of 35.5° sun incidence and 10° top spin. 
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111. Results and Discussion 

The comparison between the L’Garde STK and the S5 simulations using the percent error calculations as above 
resulted in the two simulations matching to a high degree of accuracy. Three different test cases were used to assess 
the performance of the trajectory analysis. An additional test measured the effect on performance of the maximum 
allowable attitude error, and a final comparison between the flat plate and shaped sail models indicates an impact of 
sail shape on trajectory. 

A. Flat Plate Comparisons 

The first test case was propagation using both sets of software with no SRP forces acting on the spacecraft in 
either simulation for the maneuver sequence length of thirty days. This was done to define the mismatch between the 
models in terms of the effects of differences in simulation such as integration, propagation, and gravity/force field 
models without the sail. The results of this comparison are shown in Table 1 in the first row for both methods 1 and 
2 (“Spacecraft (No SRP)”). These percent error values are in the range of approximately 0.1% error on both position 
and velocity across error calculation methods. This shows decent agreement in general between the two propagators, 
and while the accuracy could be improved, it is sufficient for the purposes of validating the L’Garde model at this 
stage of the mission design process. 

Table 1. Results of the comparison between STK and S5 in percent error on position and velocity. 


(Method 1) 

Position Percent Error 

Velocity Percent Error 

Spacecraft (No SRP) 

0.0158% 

0.0487% 

Without Sail 

40.2812% 

195.1938% 

Sail 

0.0854% 

0.3012% 

(Method 2 Ar, Av) 

Position Percent Error 

Velocity Percent Error 

Spacecraft (No SRP) 

0.0222% 

0.0921% 

Without Sail 

63.5495% 

64.0808% 

Sail 

0.1348% 

0.0989% 



5 











Figure 3. Images from STK showing the final STK (turquoise) and S5 (yellow) trajectories being 

compared. 

Additionally, the STK trajectory as propagated with the solar sail was compared to the S5 simulation with no 
SRP (no solar sail) for the 30 day period. This gave a bound on the maximum error between the simulations, 
assessed the overall effect of the solar sail on the trajectory, and also allowed the effects of a sail deployment failure 
to be estimated in terms of the percent error values (especially to compare to the conservative 20% requirement). 
The results are given in Table 1 in rows 2 for both methods 1 and 2 (“Without Sail”). As seen in the table, the 
smallest percent error value for this case is 40% on position using the first method, significantly more than the 20% 
requirement. 

Finally, the comparison between the two simulations with SRP being modeled in both STK and S5 is given in 
Table 1 in row 3 for both methods as “Sail.” This comparison was done with the Derbes and Lichodziejewski SRP 
model 10 in S5 to match the SRP model in STK. It resulted in values of percent error on the order of ~0. 1%, with the 
largest being 0.13% error on position using method 2. Also, the trajectories from this case were plotted in STK and 
can be seen in Figure 3. 

The results of errors on the attitude of the sail are shown in Table 2. These are the effects on the percent error 
values for the navigation mission requirement, and are all relative to 0% error at the trim conditions. Here, it can be 
seen that variance in the sun incidence angle had a much greater effect on navigation error than variation of the top 
spin angle. This is expected as changing the sun incidence angle affects a larger force component than changing the 
top spin angle because of the trim orientation of the sail. 

Table 2. Results of percent error from ±2° variation of top spin and sun incidence angles. 


Sun Incidence Angle 

Position Percent Error 

Velocity Percent Error 

33.5° 

4.2162% 

4.5155% 

37.5° 

4.1894% 

4.4581% 

Top Spin Angle 

Position Percent Error 

Velocity Percent Error 

8° 

0.9980% 

0.8951% 

12° 

0.9980% 

0.8952% 


B. Shaped Model Results 

As stated above, the STK model tested in the validation process was a flat plate sail using the L ’Garde optics 
model. However as can be seen in Figure lb, the L’Garde sail has a very non-flat shape. In order to model the sail 
shape, L’Garde has developed force coefficients that are similar to lift coefficients for an airplane that consist of 
non-dimensionalzed force coefficients indexed against sun incidence angle. Previously L’Garde has delivered two 
sets of these coefficients to NASA. During the In Space Solar Sail Project 14 , L’Garde delivered coefficients for a 
10,000 square meter sail. Another set of coefficients was delivered for the ST9 proposal of 2006, which were 
developed for a 950 square meter sail. L’Garde is basing the Solar Sail Demonstration project on the ST9 
coefficients, since those are closer to the 1200 square meters sail size, and also because both designs include 
spreader bars (while the In Space sail design did) and thus the ST9 sail shape is probably also more similar to the 
test sail being manufactured for the Solar Sail Demonstration project. 

In order to model the shape of the L’Garde sail, we use the Generalized Solar Model (GSM) developed for S5 15 . 
The GSSM uses tensors to model non-flat sails. In the case of the L’Garde sail, since L’Garde has produced and 
delivered non-dimensionalzed coefficients, Rios-Reyes developed a method of fitting the L’Garde coefficients to the 
GSSM tensors using a least squares fit 16 . Rios-Reyes converted the coefficients for both the ST9 and In Space 
L’Garde models. We analyzed both, but here we only report results from the converted ST9 coefficients. 
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Figure 4 presents the results of the least squares fit to the coefficients. The method developed by Rios-Reyes of 
fitting the tensors to the force coefficients does an excellent job of matching the L’Garde force coefficients across 
the entire range of sun incidence angles of interest. More details of how we modeled the L’Garde shape appear in 
the Appendix. 

The results of the S5 propagated orbit using the ST9 coefficients are presented in Table 3. These results are 
compared to the same trajectory that appears in Table 1, which is the STK flat plate Derbes and Lichodziejewski 
SRP model. The larger deviation from the STK flat plate model with the L’Garde shape model is to be expected due 
to the large deviation from a flat plate of the L’Garde design. Based on the fidelity of the model in the L’Garde 
force coefficient model, this is a fairly significant result and is published here for the first time. 



Figure 4. S5 Tensors Least Squares Fit to L’Garde Force Coefficients 


Table 3. Comparison of S5 L’Garde GSM Tensor Model to L’Garde STK Flat Plate 


(Method 1) 

Position Percent Error 

Velocity Percent Error 

Shape to Flat Plate 

3.5327% 

16.5435% 

(Method 2 Ar, Av) 

Position Percent Error 

Velocity Percent Error 

Shape to Flat Plate 

5.5657% 

5.4311% 


Currently a shape model for STK of the L’Garde sail based on the ST9 coefficients is under development. 
However, those results were not available at publication time and so a comparison of the S5 shape model to the STK 
model remains future work. 


V. Conclusions 

The focus of this work was to validate the L’Garde STK solar sail simulation and assess the effect of the non-flat 
shape of the L’Garde sail on meeting mission requirements. The NASA-developed simulation S5 matched the 
L’Garde STK simulation of the Sunjammer 30-day maneuver sequence very closely for the flat plate model, with 
percent error values on the order of ~0.1% on position and velocity using both methods. The degree to which the 
simulations match can be seen by observing that these errors are not much larger than those associated with the 
comparison without the sail. This shows good agreement between the sail models, but with some differences still 
remaining. While the difference between the two software tools could be improved to be even more accurate, these 
results indicate that the L’Garde simulation is accurate enough to be used for the mission requirements. Further 
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refinements may be introduced at the Flight Readiness Review (FRR) phase of a future mission, and as the L’Garde 
simulation is updated and finalized, the existing S5 model can also be modified to continue to validate that the STK 
simulation is accurate. Additionally, the result of a minimum of 40% error in the comparison of the sail model in 
STK to the simulation with no SRP (no sail) in S5 implies that while the 20% requirement is not exceedingly 
demanding, it has been determined that this requirement cannot be met if the sail fails to deploy. 

The attitude error study shows that the 2% requirement to point the sail will not cause the position requirement 
of 20% to fail even with the maximum error. The maximum deflections from trim for each angle were studied, with 
results indicating that errors in attitude will have a significant effect on the navigation requirement, and also that as 
expected, errors in the sun incidence angle have a greater effect than errors in the top spin angle. 

The results of the shaped model comparison show that the shape of the L’Garde sail, as expected, has a large 
effect on the trajectory. This result indicates that any mission plans with a L’Garde sail, or a sail having a similar 
large deviation in shape from a flat plate, should incorporate a shape model for accuracy. The results also imply that 
calibration of a non-flat sail early in a future mission will be important, since the differences between the model and 
actual thrust are likely to be greater for a sail that deviates significantly from a flat plate. 

The second error method is recommended for use when comparing the L’Garde simulation to the as-flown 
trajectory, as it resulted in values of percent error that are more consistent between velocity and position. This 
recommendation has been passed to the Solar Sail Demonstration Project office. A more complex analysis of the 
orbit propagation errors could have been performed, but the project requirement was written as a percent deviation 
in the final position of the sailcraft, so the analysis was performed on that basis. 

Future work could extend the analysis to a comparison of the shape model being developed for STK and also other 
shaped models that publish results. For the attitude errors, a statistical investigation of the effects of combined error 
(varying both top spin and sun incidence together) on the attitude of the sail could provide valuable information for 
use during the comparison of the L’Garde simulation to the as-flown trajectory. A more complex analysis of the 
trajectory errors could also be performed, perhaps incorporating a typical covariance matrix for a notional 
navigation filter. A significant finding implied by the shape model results is that even flat solar sails need to 
consider modeling things such as boom bending or areal differences between sail quadrants as a sanity check on the 
flat plate model. This work will help inform future missions for the L’Garde sail as well as other current sail 
projects such as the NASA NEA Scout and Lunar Flashlight missions, the Surrey Cubesail, or the Planetary Society 
Lightsail. 
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Appendix 

The constants from the SRJP force model developed by Derbes and Lichodziejewski.'°This is the model used by 
STK for the flat plate comparison. 

R p = 0.85268583078932 
9 p = 25.9315208633397 
p d = 0.894100805958922 
£ f = 0.05 
e b — 0.62 
B f = 2/3 
B b = 2/3 

For the modification to the sine term in Equation (3) we used: 
k= 1.037260834533588 (i.e., 9 p — kdi for the sine term) 

For the shape model in section IV, we used the following equation which is Rios-Reyes 1 6 Equation (15): 


F c = 


(i 2 t i 0 —Apsr\r 3 0 

0 a2'/'3 —2 psf\ —2psrl 


-a 3 7\r 3 

-a 3 r| 


For which the tensors are: 


4 

4 

4.3 

4.3 

A 


Jii = 0.1172 
JIjt 2.9733 
Jl 13 = 0.0042 
/! 3 3 = 1-1347 
/|= -1.1079 

The above tensors were matched to the L’Garde force coefficients using a Least Squares Fit by Rios-Reyes 16 . 
And the optical coefficients are (Mclnnes model): 


s = 1.0000 
p = 0.8941008 
£f = 0.05 
£ b = 0.40 
B f = 0.79 
B b = 0.55 

The above coefficients are what Rios-Reyes fit the tensors to, and are his approximation of the Derbes model. These 
coefficients are also reduced by Rios-Reyes to: 


a 2 = B f {l-s)p + (1 -p) 
a 3 = 1 — ps = 0.10590 


e f B f -e b B b 
. £ f~ £ b 


-0.04248 


r lt r 2 ,f 3 represent the sail to sun vector in the body frame of the sail where the 3- axis is perpendicular to the normal 
of the sail surface (see Figure 2). 
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